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ABSTRACT 

The effect of a magnetic field on the linear phase of the advective-acoustic instability is investi- 
gated, as a first step toward a magnetohydrodynamic (MHD) theory of the stationary accretion shock 
instability taking place during stellar core collapse. We study a toy model where the flow behind a 
planar stationary accretion shock is adiabatically decelerated by an external potential. Two magnetic 
field geometries are considered: parallel or perpendicular to the shock. The entropyvorticity wave, 
which is simply advected in the unmagnetized limit, separates into five different waves: the entropy 
perturbations are advected, while the vorticity can propagate along the field lines through two Alfven 
waves and two slow magnetosonic waves. The two cycles existing in the unmagnetized limit, advec- 
tiveacoustic and purely acoustic, are replaced by up to six distinct MHD cycles. The phase differences 
among the cycles play an important role in determining the total cycle efficiency and hence the growth 
rate. Oscillations in the growth rate as a function of the magnetic field strength are due to this varying 
phase shift. A vertical magnetic field hardly affects the cycle efficiency in the regime of super- Alfvenic 
accretion that is considered. In contrast, we find that a horizontal magnetic field strongly increases 
the efficiencies of the vorticity cycles that bend the field lines, resulting in a significant increase of the 
growth rate if the different cycles are in phase. These magnetic effects are significant for large-scale 
modes if the Alfven velocity is a sizable fraction of the flow velocity. 

Subject headings: instabilities - - shock waves - - magnetic fields — magnetohydrodynamics — 
supernovae: general 



1. INTRODUCTION 

The delayed energy deposition by neutrinos may 
be a viable explosion mechanism of massive stars 
(|Bethe fc Wilsonf [1985) if multidimensional hydrody- 
namical instabilities are taken into account during the 
first second after s hock bounce, in the inner 200km 
of th e massive star (jBuras et~aTT[200l iMarek fc Jankal 
2009). By capturing some gas in large scale convective 
cells, these instabilities increase the e xposure time of 
this gas to the heating by neutrinos (IMarek fc Janka 
20091: Murphy fc Burrowsll2008l; IFernandez fe Thompson 



2009a). The I = 1,2 deformation of the shock seems 
to be mainly due to the Stan ding Accretion Shoc k 
Instability (SASI), discovered bv IBlondin et al.l (|2003f ). 
This asymmetric instability may influence th e kick 
and spin of the residual pu lsar (|Scheck et al.l 120061 : 
IBlondin fc Mczzac appal 120071). tri gger the excitation 
of g-modes ((Burrows et al. 20061) and the emissio n 
of gravitational waves (|Marek et afl 120091 : lOtH 120091 ) . 
These potential consequences of SASI are based on 
numerical simulations where the magnetic field is 
neglected, because the magnetic pressure is thought 
to be small compared to the thermal pressure in 
a typical core collapse. Other models of core col- 
lapse supernovae, in which magnetic forces play a 
dynamical role, usually rely on a fast rotation rate 
from which the magnetic energy is extracted by 
field w inding and loca l insta bi lities (|Akivama et alj 
120(31 iMoiseenko et~aTl l200l iShibata et all 120061: 
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iBurrows et alj 120071: lObergaulinger et al.l 120091 ). Sur- 
prisinglv. lEndeve et al. I (120081) argued that the turbulent 
flow associated with SASI may be able to amplify an 
initially weak magnetic field to a dynamically significant 
strength of 10 15 G, even if the progenitor is not rotating. 
This latter result, if confirmed, shows the importance of 
understanding SASI in the context of MHD even if the 
progenitor is weakly magnetized. 

As a first step in this direction, we investigate the 
linear phase of SASI in a moderate magnetic field. 
The linear mechanism of SASI has been iden ti fied a s 
an advective-acoust i c cycle bv IBlondin et al.l (T2 003) ; 
lOhnishi et alj (l2006l):lFoglizzo et alj (120071) :lScheck et alj 
( 2008): IFernandez fc Thompson! (|2009bf )7 The interplay 
of acoustic waves and advected perturbations of en- 
tropy /vortk2itv_lias_been illu strated using a s imple toy 
model bv lFoglizzol (|2009f ) and lSato et all (|2009f ). in which 
the flow is planar and adiabatic. A similar cycle is ex- 
pected to exist in a MHD flow, but should be modi- 
fied by the magnetic field depending on its strength and 
direction, and could potentially involve entropy waves, 
fast and slow magnetosonic waves, and Alfven waves. 
The local stability of an isolated fast MHD s hock has 
been established by iGardner fc Kruskall (|1964f ). and its 
interaction with MHD waves has been characterized by 
iMcKenzie fc WestphaH (|1970l ). However, the global sta- 
bility of the cycles resulting from the linear coupling of 
these MHD waves through a region of deceleration be- 
low the shock has never be en studied. For this purpose 
we extend the toy model of [Foglizzo ( 2009) to a magne- 
tized fluid, and investigate two possible orientations of 
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the magnetic field. 

It should be noted that the existence of SASI type 
instabilit ies in a magnetized flow has been recently pro- 
posed bv lCamus et al.1 ()2009f ) to explain the oscillations 
of the termination shock in their MHD simulations of the 
relativistic pulsar wind in the Crab nebula. The present 
study can thus be considered more generally as a first 
step toward understanding the MHD extension of the 
classical advective-acoustic cycle. 

In Sect. 2 we describe the stationary flow of our MHD 
toy model and define the method of our linear analy- 
sis. The efficiencies of all existing cycles are calculated 
in Sect. 3, and the growth rate of the resulting eigen- 
mode is computed in Sect. 4. The results are summarized 
and discussed in Sect. 5. Most algebraic derivations have 
been separated from the main text for the sake of clarity, 
and can be found in the Appendices. 



2. SET UP: THE MAGNETIC TOY MODEL 

2.1. The stationary flow 

In the toy model of the adv ective-acoustic instability 
introduced by Foglizzo (|2009f ). an ideal gas character- 
ized by an adiabatic index 7 = 4/3 flows in the nega- 
tive z-direction. This planar adiabatic flow is decelerated 
through a stationary shock at z — z sn , with an incident 
Mach number A^i. It is further decelerated over a dis- 
tance 7?v < H in the vicinity of zv = -z s h — H, by a 
decrease A$ of the external potential : 



tanh &7f) +1 



(i) 



We add to this flow a magnetic field (B x ,0,B z ) that is 
oriented either vertically (B x = 0), or horizontally (B z = 
0), as illustrated in Fig. [TJ This magnetic field does not 
depend on x or y. 

Let us denote by P the gas pressure, p the density, 
c = \f~(Pjp the sound speed, and v the fluid velocity. 
S = [log ((P/P sh ) / (p/psh) 7 )] / (7 - 1) is a dimensionless 
measure of the entropy, vm = \J Bf / pop is the Alfven 
speed along the direction i = x or z . 

The stationary flow is governed by the conservation of 
the mass flux pv z , the energy density v z /2 + c 2 / (7 — 1) + 
v Ax + an d the entropy S. The conservation of B z and 
Bx/ P follows respectively from the conservation of the 
magnetic flux and the frozen-in field lines in ideal MHD. 

A generalization of the Rankine-Hugoniot jump con- 
ditions to MHD characterizes the post-s hock flo w as a 
function of the upstream conditions (Eqs. (|A1IIA2[) in Ap- 
pendix A). The additional pressure due to a horizontal 
magnetic field decreases the compression at the shock 
and increases the post shock Mach number by less than 

10% if UAsh < «sh- 

The flow below the adiabatic region of deceleration is 
characterized by the size A$ of the potential jump, or 
equivalently by the ratio of sound speed c- m /c ou t, using 
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Fig. 1. — Schematic view of the toy model. A magnetic field that 
is either vertical o r horiz ontal is added to the hydrodynamical toy 
model of Foglizzo (2009). Six waves propagate from the shock to 
the inhomogeneous region near (an entropy wave S, two Alfven 
waves a±, two slow waves Szfc and a fast wave f+) and couple to the 
upward propagating fast wave f— . The linear coupling between a 
wave i and the wave f— is described by the coupling coefficients 
and defined in Sec. 2.2 
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Unless otherwise specified, we use the values of the pa- 
rameters cf n /cl ut = 0.75, H^/H = 10~ 3 , Mi = 5 when 
studying the effect of the magnetic field. This choice 
implies an increase of the height of the potential jump 
with the field strength, which becomes significant for 
fAxsh/^sh > 1 — 1.5, i.e. when the magnetic pressure be- 
comes comparable with the thermal pressure (Fig. [2]). Al- 
ternatively, one could choose to keep the potential jump 
constant. This choice of parameterization is unimportant 
in the present study, which is focussed on the regime of 
weak magnetic pressure. 

2.2. Linear perturbations 
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Fig. 2. — Effect of the magnetic field on the stationary flow. 
Upper panel: Ratio of magnetic pressure to thermal pressure just 
below the shock (full line) and below the potential jump (dashed 
line). Bottom panel: Change of the potential jump if the compres- 
sion is kept constant with c? n /c^ ut = 0.75. The potential jump is 
increased by ~ 20% at i>a = 1. 



The 1-D stationary flow is perturbed in 3-D, using 
Fourier transforms in time and in the directions x and y 
where the boundary conditions are periodic. The struc- 
ture of the perturbation of a physical quantity / is thus: 
Sf(z) x expi (k x x + k v y — cot). For a given horizontal 
size of the computational domain L x and L y (we chose 
L x — L y = 4ff), the horizontal wavenumbers take dis- 
crete values ki = 2irrii/Li, where i = x,y and n, is the 
number of horizontal wavelengths in the i-direction. We 
impose a leaking boundary condition below the potential 
jump (i.e. no wave propagating upward from below), and 
leave the supersonic flow unperturbed ahead of the shock. 
The corresponding equations are detailed in Appendices 
A and B. 

In each uniform section of the flow, perturbations can 
be decomposed into 7 waves: two fast magnetosonic, two 
slow magnetosonic, two Alfven and one entropy wave. 
If the shock is a fast MHD shock, only the fast mag- 
netosonic wave is able to reach it from below. We do 
not consider a slow MHD shock because it would re- 
quire Alfven waves to be able to propagate upstream 
of the shock (which would correspond to an unrealisti- 
cally strong magnetic field), nor transAlfvenic disconti- 
nuities, which ha ve been shown to be non evolutionary 
(|Svrovatskii|[T9 59). As a consequence, we always consider 
a situation where, due to advection, 6 of the 7 waves are 
propagating downward in the post-shock region (Fig. [1]). 

A global mode can be decomposed into 6 cycles, asso- 
ciated with each of the six waves propagating downward, 
while the upward propagation always corresponds to a 
fast magnetosonic wave. Waves are linearly coupled at 
the shock and in the regi on of flow gradients around Zv • 
Following iFoglizzol (|2009D . the coupling efficiency Qy de- 



scribes the ratio of amplitudes of a perturbation SA of 
a physical quantity A which is propagated from z s ^ to 
zv, generates a fast magnetosonic wave in the region of 
gradients, which propagates back to the shock: 
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where the symbol + (— ) is used for waves propagating 
downward (upward) in the frame advected with the fluid. 
The superscript s± stands for slow magnetosonic waves, 
f± for the fast ones, a± for the Alfven waves and S for 
the entropy wave. 

The coupling efficiency at the shock Q l sh measures the 
creation of a wave i when the fast magnetosonic wave f — 
perturbs the shock: 
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The choice of the physical quantity A is driven by the 
simplicity of the resulting equations (see Appendices A 
and B). The total cycle efficiency Q\ defined as follows, 
does not depend on this choice: 



Q % = QvQlh- 



(8) 



We define Qtot as the sum of all the cycles constants, 
and Qvort = Q s+ + Q s ~ + Q a+ + Q a ~ is restricted to the 
cycles involving vorticity. 

It is instructive to study how these six cycles relate to 
the two hydrodynamical cycles, by considering the weak 
field limit. Since the acoustic cycle straightforwardly be- 
comes the fast magnetosonic cycle, its reflection efficiency 
is noted R instead of Q l+ . The advective-acoustic cycle 
splits between an entropy-fast cycle where the entropy 
wave is still advected, and four Alfven-fast, slow-fast cy- 
cles that contain the vorticity, which is no longer pas- 
sively advected. In the setup we consider, the entropy 
cycle typically has an efficiency of Qs ~ 1 while the vor- 
ticity cycle has an efficiency of Q V ort ~ 0.3 in the hydro- 
dynamical limit. How is vorticity distributed between 
the Alfven and slow waves ? The distinction between 
slow and Alfven waves at the shock can be anticipated 
by remembering that the velocity and magnetic field per- 
turbations of Alfven waves are along the direction k x B, 
while that of slow waves is in the plane (k, B) (or equiv- 
alently the vorticity of slow waves is along k x B, while 
that of Alfven waves is in the plane (k, B)). This ques- 
tion is discussed in the next two subsections, depending 
on the orientation of the magnetic field. In the weak 
field limit, the efficiencies associated with Alfven and 
slow waves are independent of the direction of propa- 
gation (Q a+ ~ Q a ~ and Q s+ ~ Q s ~), because both the 
Alfven and slow waves are then simply advected and the 
(±) components are undistinguishable. 

2.3. Vertical magnetic field 

A vertical magnetic field does not influence the station- 
ary flow as the fluid flows along the field lines without ex- 
periencing any magnetic force. It does however influence 
the evolution of perturbations that involve some trans- 
verse motion, in particular enabling vorticity to propa- 
gate. The equations governing the perturbations are de- 
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scribed in the Appendix B. With the magnetic field in the 
z-direction, the directions x and y are equivalent. With 
no loss of generality we choose the x-direction parallel 
to the transverse wavenumber (k y — 0) and solve a 2D 
problem in the plane (x, z). The boundary conditions 
at the shock (Eqs (B13-B14)) impose 5v y — 5B y = 0. 
The only waves along this direction would be the Alfven 
waves, which displacement is along k x B (i.e. along y 
with our choice of axis) . As a consequence, the shock os- 
cillations cannot create Alfven waves, and there are only 
4 cycles in this particular field geometry. Note that the 
absence of Alfven waves is not a consequence of choosing 
k y = 0. 

In a stationary accretion flow in uniform vertical mag- 
netic field, the ratio of the flow velocity to the Alfven 
velocity decreases downward, scaling like v/va oc p -1 / 2 . 
The transition from a superAlfvenic flow to a subAlfvenic 
flow is named the Alfven surface. Perturbations cannot 
be treated in the framework of ideal MHD in the vicin- 
ity of this surface because their wavelength becomes in- 
finitely small . Alfven waves can accumulate and be am- 
plified there (|Williamsl ll975). For the sake of simplicity, 
we restrict the present study to weak enough magnetic 
fields so that the flow is superAlfvenic everywhere. Given 
our choice of parameters, the superAlfvenic condition re- 
quires WAzsh/wsh < 0.68. The study of transAlfvenic flows 
is the subject of a separate study (Guilet et al., in prep.). 

2.4. Horizontal magnetic field 

The equations governing the perturbations in a hori- 
zontal magnetic field are described in the Appendix A. 

In general, both Alfven and slow waves are created at 
the shock. They propagate mainly along the horizontal 
magnetic field lines (slow waves also propagate slightly 
along the perpendicular direction). 

The two particular cases k x — and k y — are simpler 
because the evolution of perturbations is then planar: 

(i) if k y = 0, the boundary conditions at the shock 
impose 6v y = 5B y = while v x and B x are perturbed 
by the shock oscillations. As k x B is the y direction, 
the shock oscillations do not create Alfven waves but do 
create slow waves. 

(ii) if k x — (k_LB), the magnetic field lines are not 
bent and the only magnetic force is the magnetic pressure 
which adds up to the thermal pressure in the plane {y, z) 
perpendicular to the field. The evolution of perturbation 
is thus similar to the hydrodynamical limit, with a mod- 
ified pressure. Vorticity is advected by the flow together 
with entropy. This vorticity wave should be viewed as 
an Alfven wave rather than a slow wave, because it is in 
the direction of the magnetic field. 

3. THE COUPLING EFFICIENCY IN THE PRESENCE OF A 
MAGNETIC FIELD 

In this section we study the effect of the magnetic field 
on the coupling efficiencies, which are computed by in- 
tegrating numerically the differential system governing 
the perturbations, and using the boundary conditions 
given in the appendices. We consider plane waves with 
a real frequency which coincide with the real part of the 
eigenfrequency of a given eigenmode. Growth rates are 
discussed in Sect. 4. 

3.1. Vertical magnetic field 
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Fig. 3. — Upper panel: Efficiency of the slow cycles Q B ± as a 
function of the height of the potential jump H^j in the presence 
of a vertical magnetic field. The grey line represents the efficiency 
in the absence of a magnetic field, the black lines correspond to 
v Azsh = 0.5t) s h. The dashed line is the (-) wave, the full line the 
(+) wave. The wave propagating up (Q B — ) is more affected by 
the cutoff than the wave propagating down (Q s +). Bottom panel: 
The cycle efficiencies are normalized by Qo defined as the limit of 
this efficiency when the potential jump is compact. The cutt off 
happens when k z H\j > 1 — 2, where fc z is the vertical wavenumber 
of the slow wave involved in the cycle. 

By solving the evolution of perturbations described in 
Appendix A, we remark that the influence of a vertical 
magnetic field on the individual coupling constants of the 
different cycles is moderate (< 20%) in the long wave- 
length limit, where the compact approximation is valid. 
The efficiencies of the two slow cycles at shorter wave- 
length differ in a way whic h can be under stood in terms 
of the cutoff described by iFoglizzol ([2009) . In the weak 
field limit, the vertical wavenumber of slow waves can be 
approximated by: 

kf ~ -J^—. (9) 

V z T v Az 

The slow wave propagating against the stream (-) has a 
shorter wavelength than the wave propagating with the 
stream (+). It is thus more sensitive to the cutoff induced 
by the finite scaleheight of the gradients (Fig. [3l upper 
panel). As in the hydrodynamic case, the cutoff takes 
place for k z H\j ~ 1 (Fig. [3J bottom panel). 

3.2. Horizontal magnetic field 
3.2.1. Straight field lines (k _L BJ 

We observe the presence of two slow cycles, which 
were absent in the non- magnetic case (Fig. [4j. These 
slow waves taken together (+ and - are equivalent here 
as these slow waves do not propagate) do not contain 
any velocity or vorticity perturbation: they are simply 
a perturbation of the magnetic field strength and den- 
sity in such a proportion that the total pressure is un- 
perturbed. As their efficiency depends weakly on the 
transverse wavenumber, we focus on the simplest case 
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Fig. 4. — Apparition of a slow cycle in the presence of a horizontal 
magnetic field and k_LB. This cycle carries no vorticity. It is only a 
perturbation of the magnetic field and density in such a proportion 
that the total pressure is not perturbed. 



of n y = 0. In this case, the presence of the slow waves 
can be explained using the conservation of magnetic flux, 
which perturbation 5A = 8B X / B — dp/ p should vanish 
at the shock (Eqs. (|A"29|| ). 

While fast and Alfven waves do not perturb the mag- 
netic flux, the entropy wave is responsible for a perturba- 
tion 5 As = (7 — 1) SS/j. Slow waves have to be created 
at the shock in order to keep 6A s h = 0. 

In the deceleration region, the slow waves necessarily 
disturb the pressure equilibrium and create an acous- 
tic feedback, in order to satisfy the conservation of both 
the mass flux and the magnetic flux across the poten- 
tial jump. This acoustic generation is very similar to the 
acoustic feedback produced by the deceleration of the en- 
tropy wave, required by the conservation of entropy and 
mass flux across the potential jump. 

The effect of the magnetic field on other cycles is only 
minor and depends on the parameterization of the po- 
tential jump. For example the efficiency of the entropy 
cycle is slightly increased if the compression is kept con- 
stant when the magnetic field is varied, while it is slightly 
decreased if the height of the potential jump is kept con- 
stant. 

3.2.2. Bending the field lines: amplification of the vortical 

cycles 

If k || B, the slow cycles carrying the vorticity are 
strongly amplified by the presence of a horizontal mag- 
netic field (Fig. [SI upper panel). The amplification is al- 
most linear with the magnetic field strength and reaches 
a factor ~ 5 at VAxsh = Vsb- Note that this effect is more 
important than the one described in Sect. 3.2.1, which 
efficiency only reached Q s± ~ 0.1 at VAxsh — w s h- 

When k is oblique with respect to B, the vorticity is 
distributed in both the Alfven and the slow waves. The 
Alfven cycles are also amplified, although to a lesser ex- 
tent than the slow cycles (Fig. [5j bottom panel). When 
varying the angle between k and B, the total vortical cy- 
cle efficiency (the sum of the four slow and Alfven cycles) 
increases with k z /k. 

Although we do not provide an analytical description 
of this amplification, one can gain insights into its nature 
by first remarking that this effect is not restricted to slow 
waves and hence is not due to the compressional nature of 
the slow waves. Second, this amplification appears only 




Fig. 5. — Amplification of the vorticity cycles in the presence of 
a horizontal magnetic field. Upper panel: k || B (n x = 1, L x = 4), 
only slow waves are created at the shock. Bottom panel: k oblique 
with respect to B (n x = n y = 1, and L x = L y = 4v2). Both slow 
and Alfven waves are created at the shock, and amplified by the 
presence of the horizontal magnetic field. 



when the field lines are bent by the perturbations. One 
should however note that a vertical magnetic field bent 
by perturbations does not lead to the same amplification. 

4. BUILDING OF A GLOBAL MODE BY ADDING UP 
CYCLES 

4.1. Method 

The cycle efficiencies can be calculated using a com- 
plex frequency to = uj r + icJi, where to r is the oscillation 
freque ncy and ujj the g rowth rate. Generalizing the anal- 
ysis of iFoglizzol (j2009f ) . the complex eigenfrequency u of 
a global mode obeys the following relation: 



= Q s + Q s+ + Q s - + Q a 



R = l. (10) 



The cycle efficiencies computed in Sect. 3 with a real 
frequency can be used to characterize marginal stability 
(IQtotl = 1)- |Qtot I < 1 indicates a stable mode, and 
|Qtot| > 1 an unstable one. 

An estimate of the growth rate associated wit h a single 
cycle in the limit uji <C uj r has been given in iFoglizzo 

bg|Q(w r )|, (11) 

^cyclc 

where r cyc i c is the cycle timescale. In the case of an 
advective- acoustic cycle, 



Tcyclc 



a* + M S h 

; Ai(l + M sh ) 



(12) 
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where \i is defined by [i 2 = 1 — k 2 c 2 (l — M 2 ) /uj 2 and 
Taac is the radial advective-acoustic time: 

r aac = i//(|t; sh |(l-X sh )). (13) 

The situation is more complicated if several cycles with 
different timescales are involved in the instability. One 
may however propose an approximate relation: 

~ — log |Qtot(w r )|, (14) 
T e fi 

where r c g is an effective timescale that would be an av- 
erage of the different cycle timescales. In the next two 
subsections we use Eq. f|12|) as a proxy for T e g and show 
that this is a good approximation. 

4.2. Vertical magnetic field 

Fig. [6] shows the evolution of the most unstable n x = 1 
and n x = 8 modes with a vertical magnetic field. The 
estimate of the growth rate using Q tot (Eq. (|14p) is in 
excellent agreement with the eigenvalue if n x = 8, while 
the agreement is only reasonable if n x = 1 (in which case 
the condition <C 0J r is less justified). This justifies 
a posteriori our estimate of the effective timescale T e s , 
and indicates that it is weakly affected by the magnetic 
field. Some effect of the magnetic field could have been 
expected since the vorticity is able to propagate through 
slow waves, but the (s-) wave is decelerated while the 
(s+) wave is accelerated, so that both effects on r c g can- 
cel each other to first-order. Furthermore, the most effi- 
cient cycle in our toy model is the entropic-acoustic cycle 
and at the magnetic field strength considered, the prop- 
agation speed of the fast magnetosonic waves is close to 
the acoustic one. 

The growth rate as a function of the magnetic field 
strength shows oscillations (Fig. [6j upper panel), which 
are reminiscent of the oscillations in the SASI eigenspec- 
trum observed bv iFoglizzo et al.1 (|2007T ) (Fig. 7). These 
were interpreted as the consequence of a purely acous- 
tic cycle interacting either constructively or destructively 
with the advective-acoustic cycle. Similarly the present 
oscillations can be interpreted as the constructive or de- 
structive effect of the vortical cycles (the two slow cycles). 
The phase of their cycle efficiencies Q s+ and Q s ~ varies 
with the vertical wavevector of the slow waves as k^H 
(neglecting a variation of the phase shift during the cou- 
pling process). Equation © shows that the phases of 
Q s+ and Q s ~ vary in opposite directions. This results in 
oscillations in the value of Q VO rt as well as Qtot (Fig. [6]). 
As wa approaches v zi the phase of Q s ~ varies faster and 
faster, which can recognized in the faster and faster os- 
cillations in Fig. |51 The slower phase variation of Q s+ is 
responsible for the slower background oscillations in the 
growth rate of the n x = 8 mode. 

Using a first-order expansion of Eq. ([9]), an estimate 
of the first minimum of the growth rate can be deduced 
from the criterion (k z — k z o) H = n: 

VAzsh nv sh 1 , , 

v sh ~ ujH ~ 2h (1 - M sh ) ' [ ' 

where the second estimate is obtained by approximat- 
ing the h-th harmonics frequency by ui = 2ft,7r/r aac . 
This gives VAzsh/v s h ~ 0.83 for the mode n x = 1, and 
"Azsh/Vsh ~ °- 10 for n * = 8 - 
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Fig. 6. — Effect of a vertical magnetic field on the modes % = 
1 (black) and n x = 8 (grey). Upper panel: Growth rate (full 
line) as a function of the vertical magnetic field strength. The 
quantity log (Qtot) /t is overplotted with dotted lines. Bottom 
panel: Efficiency of the vorticity cycles Q VO r- The oscillations are 
due to interferences between the downward propagating (+) and 
upward propagating (-) slow waves. 



Because the entropy and the vorticity cycles are 
roughly in phase in the non- magnetic limit, the phase 
variations due to the magnetic field induce an overall de- 
crease of Qtot and (Figs.[H]and[7]). For any strength of 
the magnetic field, some modes benefit from a construc- 
tive interference of the different cycles and are as unsta- 
ble as in the absence of a magnetic field (e.g. n x — 5, 10 
in Fig. [71). This phase effect thus never stabilizes com- 
pletely the compact toy model {H\j — 0), but is responsi- 
ble for an irregular eigenspectrum and lowers the average 
growth rate. If the deceleration region is sufficiently ex- 
tended vertically, the cutoff effect verified in Sect. 3.1 can 
stabilize high frequency modes. A strong enough mag- 
netic field may thus be able to completely stabilize the 
flow if the frequency of the constructive cycles exceeds 
the frequency cutoff induced by the finite size of the de- 
celeration region. 

4.3. Horizontal magnetic field 

Similarly to the case of a vertical magnetic field, the 
growth rate is well described by Eq. (fH| (Fig. 0). The 
oscillations are due to the propagation of vorticity along 
magnetic field lines, which changes the vertical structure 
of the slow and Alfven waves as follows: 

kf * kt = ± ^ (16) 
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Fig. 7. — Effect of a vertical magnetic field on the eigenspectrum 
of the compact toy model (Hy = 0). The most unstable mode 
is shown for each n x from 1 to 10. White diamonds are without 
magnetic field, black diamonds with ^Azsh = 0.5u s h. 

Transverse slow modes are most affected by the magnetic 
field. The magnetic field strength of the first minimum in 
the growth rate oscillations, due to the destructive effect 
of the slow or Alfven cycles on the entropic cycle, can be 
estimated as: 



L, 



k x H 2n x H 



(17) 



This gives UAxsh/vsh ~ 2 if n x = 1, and VA xs h/v s h ~ 0.25 
if n x = 8. 

Furthermore, due to the amplification of the vortical 
cycles described in Sec. 13.2.21 the growth rate at the 
oscillation maximum increases with the magnetic field 
strength (Fig. [5]). As an example, the mode n x = 8 
grows four times faster at i?Axsh — v s h than at VAxsh = 0. 
The eigenspectrum is more irregular than the hydrody- 
namical one, due to the interferences between the dif- 
ferent cycles, and shows higher maximum growth rates: 
WjTaac ~ 0.8 if VAxsh = «sh, compared to WiT aac ~ 0.4 if 
VAxsh = (Fig. El). 

5. DISCUSSION AND CONCLUSION 

The main results of this study can be summarized as 
follows: 

1. In the presence of a magnetic field, the advective- 
acoustic cycle splits in up to 5 different cycles. The 
acoustic wave becomes a fast magnetosonic wave, 
while the entropy- vorticity splits into an entropy 
wave, 2 Alfven waves and 2 slow magnetosonic 
waves. 

2. The acoustic cycle becomes a fast magnetosonic cy- 
cle. Its efficiency is not significantly affected by the 
presence of a magnetic field in any of the configu- 
rations considered in this paper. 

3. The propagation of the vorticity through slow and 
Alfven waves leads to a phase difference between 
the different cycles, which interact either construc- 
tively or destructively depending on the mode con- 
sidered. The consequence of these interferences is 
a more irregular eigenspectrum. 




BA v /v c ' 6 

Ash Bh 



Fig. 8. — Same as Fig. [6] but with a horizontal magnetic field. 
The modes presented here have k || B: n x = 1 — 8, n y = 0. 




Fig. 9. — Effect of a horizontal magnetic field on the eigenspec- 
trum. The most unstable mode is shown for n x varying from 1 to 
10, and n y = 0. White diamonds are without magnetic field, black 
diamonds with t'Axsh = v sh- 



4. In the superAlfvenic regime that we investigate, a 
vertical magnetic field hardly changes the coupling 
efficiencies in the compact toy model. It only af- 
fects the cutoff associated with the size Hy of the 
deceleration region by changing the vertical struc- 
ture of the slow waves. 

5. The vortical cycles are strongly amplified by a hor- 
izontal magnetic field if the field lines are bent 
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(k.B 7^ 0). Both the Alfven cycles and the slow 
magnetosonic cycles are amplified. This amplifi- 
cation of the vortical cycles leads to faster growth 
rates of the modes with n x ^ 0. 



These results confirm that the mechanism of the 
advective-acoustic cycle can be generalized to a magne- 
tized environment, such as core collapse supernovae or 
possibly the termination shock of pulsar wind nebulae. 
In order to estimate the significance of these results for 
SASI during core collapse, we can estimate the magnetic 
field strength needed to affect the advective-acoustic cy- 
cle. The desynchronization effect studied in Sects. 14.21 
and 14.31 requires «Ash ~ ^sh for low frequency, low n x 
modes. Eq. (fT7|) can be extrapolated to a spherical core 
collapse as follows: 



conservation of the magnetic flux 



VA 

v 



2y/l(l + \)H 



(18) 



For the dominant 1 = 1 — 2 modes, with i? s h = 110km 
and H = 50km, the desynchronisation of the cycles is 
expected for va/v ~ 0.4 — 0.8 . 

The magnetic field strength at which the amplifica- 
tion of the vortical cycles becomes significant is typically 
VAsh ~ (0.5 — 1) x i; sh in our toy model. A more detailed 
study should determine whether this amplification takes 
place at the shock or in the region of deceleration. 

Interestingly the magnetic effects described in this pa- 
per seem to be significant when the Alfven speed is com- 
parable with the advection speed rather than when the 
magnetic pressure is comparable with the thermal pres- 
sure. These two criteria are related through the Mach 
number M.: 



mag 
Pth 



2 V v J 



(19) 



In the subsonic flow the pressure ratio is thus smaller 
than va/v, by a factor ~ 0.06 at the shock if A^sh ~ 0.3, 
and as low as 10~ 3 if ~ 0.05 at the coupling radius. 
This raises the possibility that the magnetic field may 
affect the growth of SASI through the magnetic tension, 
even if the field is so weak that the magnetic pressure is 
negligible. 

The strength of the magnetic field present in the iron 
core of a star before the collapse is very uncertain. The 
best estimate to d ate is: B$ ~ 5.10 9 G and B r ~ 10 6 G 
(|Heger et a l. 2005), which is too weak to have any effect 
on SASI: indeed, flux conservation during the collapse 
(B/p 2 / 3 = cste, e.g. iShibata et ail (120061 ) ) would lead 
to WAsh/wsh ~ 0.001 at the shock (r s h ~ 150km, p s h ~ 
10 9 g.cm -3 ), and va/v ~ 0.025 at the proto-neutron star 
surface (rpNs ~ 50km, ppns ~ 10 n g.cm~ 3 ). 

A reference Alfven speed near the surface of a strongly 
magnetized proto-neutron star can be estimated by as- 
suming that the fossil magnetic field is strong enough 
to give birth to a magnetar (B ~ 5.10 15 G) by simple 



^-0.93 x 



B 



NS 



5.10 15 G 



4.10 14 g.cm- 



3\ 2 /3 



7/6 



Pns 
0.3MR.S" 1 



(50km) 



Such a ffetfcc&fa^e 3 enough tcMffect SASI significantly 
if the amplification of the vortical cycle takes place at 
the coupling radius, which is slightly a bove the proto- 
neutron star surface (jScheck et al.l 120 08). Conversely, if 
the amplification takes place at the shock, even a fossil 
magnetic field corresponding to magnetar strength would 
have a negligible effect on SASI (wAsh/wsh ~ 0.027). Note 
however that even stronger magnetic fields are sometimes 
considered in the lite rature, with Alfven waves propagat- 
ing above the shock (jSuzuki et al.ll2008h . 

Given the simplicity of our toy model, we cannot ex- 
pect it to capture more than the first-order effects of the 
magnetic field on the advective-acoustic cycle. More ac- 
curate estimates should consider a more realistic setup 
that includes radial convergence and non-adiabatic ef- 
fects, which are expected to affect the relative impor- 
tance of the entropy and vorticity cycles. These effects 
could be estimated in the c ylindrical geometry used by 
lYamasaki fc Foglizzo! (|2008ft . 

The role of the Alfven surface has not been considered 
in the present study, because the region of acoustic feed- 
back was assumed to be superAlfvenic for the sake of 
simplicity. If the magnetic field were vertical and strong 
enough, the magnetic extension of the advective-acoustic 
cycle to MHD cycles would have to take into account the 
possible coupling processes taking place at the Alfven 
surface. The physics of trans Alfvenic accretion flows will 
be considered in a forthcoming study (Guilet, Fromang 
& Foglizzo, in prep.). 

The amplification of th e magnetic field obser ved in the 
numerical simulations by lEndeve et al.l (|2008l ) seems to 
take place in the nonlinear regime of the MHD-SASI in- 
stability, and is thus beyond the scope of the present 
linear s tudy. The top o logy o f the magnetic field consid- 
ered by lEndeve et al.l (|2008[) is initially vertical (a split 
monopole), for which our analysis did not reveal any 
significant magnetic destabilization. Based on the cur- 
rent understanding of the nonli near saturation pro cess 
of SASI by parasitic instabilities (|Guilet et al.ll2009l ). we 
would rather anticipate a larger saturation amplitude of 
SASI if the field were horizontal near the shock, due to i) 
the higher growth rate of the cycles involving the bend- 
ing of field lines by vorticity perturbations, ii) the re- 
sistance of magnetic tension to the development of par- 
asitic instabilities such as Rayleigh- Taylor and Kelvin- 
Helmholtz instabilities. As noted in Guilet et al. (2009) 
however, the Rayleigh- Taylor instability may still be able 
to develop as a parasitic instability in the direction per- 
pendicular to the magnetic field, without bending the 
field lines. Besides, the presence of other magnetically 
induced parasitic instabilities cannot be ruled out. Al- 
together, these qualitative statements are insufficient to 
explain the g r owth of the magnetic energy observed by 
lEndeve et all (|2008l ). 

J.G. is thankful to the Physics Department of the Uni- 
versity of Central Florida for its hospitality. This work 
has been partially funded by the Vortexplosion project 
ANR-06-JCJC-0119. The authors are grateful to the 
anonymous referee for his comments that improved the 
clarity of the paper. 
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APPENDIX 
HORIZONTAL MAGNETIC FIELD 

Stationary flow 

Denoting by Mai — vai/c\ the ratio of the Alfven speed and the sound speed ahead of the shock, the postshock 
Mach number is influenced by a horizontal magnetic field as follows: 



M\ 



X 



1 



M 2 



Al 



-( 7 -l)A4! lX 3 -( 7 -l)M 



where x = Psh/pi is the compression factor, given by: 



(Al) 



X 



(1 + + lM 2 Ai y + (2 - 7 ) ( 7 + 1) M\M 2 A1 (1 + + iA4|0 



(2-j)M 2 



(A2) 



Al 



Definition of the variables 

The 7 variables defined in Eqs. (| A7B A9|) are chosen to describe the perturbations. This choice is guided by the quest 
for the simplest possible differential system in which the stationary flow gradients do not appear. These variables have 
the advantage of being conserved through a compact potential jump. 6 A is the perturbation of B/ p, which is conserved 
in the stationary flow due to magnetic flux conservation. Sf is the perturbation of the Bernoulli invariant and SS 
is the entropy perturbation. SE X and SE y are the perturbation of the transverse electric field along the directions x 
and y, normalized by the stationary electri c field E v = —vB. SK\ des cribes the velocity along the y direction and is 
conserved in the absence of magnetic field ([Yamasaki & Foglizzo 20081). 

5B X Sp 



6A=- 

5E x i 
Sv x i 

bE y , 
5Ki\ 

8Si 



B 

:VSV Z ~ 

_5_By 

B 

:6V X ~ 

5B X 
B 

: iuSv,, 



P 

2cSc 



7 



1 



, SB X 
' B 



Sp 
P 



v 



SB, 



7-1 



v B 

Sv z 

V 

- ikySf 
S_P_ 



5_p 
P 



7- 



(A3) 

(A4) 

(A5) 

(A6) 

(A7) 
(A8) 

(A9) 



The usual physical quantities can be expressed as a function of the new variables through the following set of 
equations: 



■SE„ 



-SE X 



= 6i 



SB X 
B 

S_By 

B 

5B Z 
B 

Sv x 

SVy 

Sv z 
v 

S_p 

P 
Sc 2 

9- = (7-l) 



Sv z 
v 



(k x SE y - k y SE x ) . 



(k x SE y 



k y SE x ) , 



: - (kySf - iSKx) , 



l + M 2 A - 
= 5E„ - 5 A 



M 2 

_ ^ 
v 



s_l 

n 2 



+ 5S + (1 + M 2 A ) 5E y + (M 2 A 



1)5 A 



SS + 5E V - SA- — 



Sw x = — 



SK i + iky 



S S 

c 2 — + v A 5A 



v A ik x 5E x 



(A10) 
(All) 
(A12) 
(A13) 
(A14) 
(A15) 
(A16) 
(A17) 
(A18) 
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where w x is the vorticity along the x direction. 



Differential system 

The differential system governing the evolution of the perturbations in a horizontal magnetic field is: 



86 A ioj 



dz v 
dSf ioj 

d6E r 



dz 
dSv x 

dz 
dSE,, 



SA 



+ %K 






Sv x - 


V 




9 Sv z 
2 — + 


2 5S 
c — 


V 


7 



(k x 6E y - k y SE x 



v 2 A SA 



v 



: l ^SE x + l ^(iSK 1 -k y 5f), 

V OJV 



i hE±[ w5 v x -{v 2 +vl) (k x 6E 
ojv v ' 



V 



kySE x ) , 



iuj - ik x 
■ — 6v x H 

V V 



? -vl) S -^ + c 2 (SA 5E y ) + (1 - 7 ) c 2 ^ 



V OJV 



38 K i iui 
— — = — 5Ki 

oz v 

dSS «w 
— — = — 8b. 

oz v 



k x v\ 



OJV 



[-k y oj5v x - (fc 2 (v 2 + v 2 A ) + oj 2 ) 8E X + k x k y (v 2 + v%) 5E y ] , 



(A19) 
(A20) 
(A21) 
(A22) 
(A23) 

(A24) 

(A25) 
(A26) 



Wave decomposition in a uniform flow 



TABLE 1 

Decomposition of the perturbations into waves: horizontal magnetic field 



fast ± 



Slow ± 



SA 
Sf 
SE X 

Sv x 

SEy 

SKi 

ss 



k z v- 



■( k l+ k l) 





Sp 






p 





1 - 



'■■ k (./ 



2 j , 2,22 
Vj^kxkz uj c —k x c 



•hi 



k z v "c~ fc x c 
fcxC k'i + k'i 



c + 



Sv x 



' V cuj c (k2 + kl) 
vjk z "t-klc 2 \ 



K c , k z wi-hic 



fej+fcj 





Sp 


iky 




p 





k x c k'i+k 2 



k z ^- 


k 2 T c 2 




k x vc k-i 






c -4( 







_ M r _|_ "I ( k l+ k io, 
2c' 2 p.ui 



1± 



Alfven ± 



Entropy 



SA 
Sf 
8E X 
Sv x 

SEy 

SK! 
SS 







ik v v 
±i 



fc z (5 It; a; 
PIP" 



l.2„,\ 5iU;b 



o v _ 



(7-1) " 
K( 7 -l)+>)f 







-ik y (yl(y-l)+c 2 )< 

SS 



In a uniform flow, the perturbations can be decomposed into 7 different waves, which are the eigenvectors of the 
differential system: an entropy wave, 2 Alfven waves an d 4 m agnetosonic waves (2 fast and 2 slow). The perturbation 
associated with each of these waves is given in Table IA.41 where oj c = oj — k z v is the frequency of a wave in the 
frame comoving with the fluid, /i 2 = 1 — k 2 c 2 (l — M 2 ) /oj 2 '', and kf = oj/c x (M =F M) / (l — -M 2 ) is the non-magnetic 
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wavenumber of the sound waves. The vertical wavevector of the magnetosonic waves is obtained by numerically solving 
the fourth order polynomial dispersion relation: 



- uj 2 k 2 (c 2 + v\) + k 2 c 2 k 2 x v\ = 0. 



(A27) 



The slow waves are then distinguished from the fast waves as follows: c 2 < u> 2 /k 2 < c 2 + v\ for the fast waves, and 
< Lo 2 /k 2 < v\ for the slow ones (va < c in the case considered here). (+) and (-) waves are distinguished by 
computing the group velocity d m / dk z along the z direction. An approximate expression of k z that is valid in the weak 
field limit is also given in T able lA~4l 

The expressions of Table IA.4I are the coefficients of a matrix which, when multiplied by the following amplitude 
vector 



r± 



6Va 

c 



Sw, r 



k 2 



a± 



,6S 



(A28) 



gives the corresponding perturbations. This matrix is inversed numerically in order to determine the amplitude of 
each wave as a function of the value of the perturbations. 

Boundary conditions 

To obtain the boundary conditions at the shock, we use the conservation of energy, momentum, and mass fluxes as 
well as the continuity of the magnetic field perpendicular to the shock, and of the electric field parallel to the shock, in 
the frame of the perturbed shock. Assuming that the upstream flow is unperturbed, this gives the following boundary 
conditions as a function of the shock displacement A£: 



SA sh -- 

Sfsh : 

o~v Xsh ■■ 



5Ey S h - 

SSsh 

1 



= 0, 

, 5f sil 



_ V2 
Vl 

V2 



ik x A( ^1 



?'2 



- 1 



J A2 



V2 



IU) . j. I Vo 

— AC 1-- 

V2 V v l 



= 0, 



— A(vi 
& 



v\ 



a$ ac 

dz c 2 



Vl 



(A29) 
(A30) 
(A31) 
(A32) 

(A33) 

(A34) 
(A35) 
(A36) 



Below the potential jump, we use a leaking boundary condition, i.e. no wave propagates upward. For this purpose, at 
z = zy — 3i?v where the flow is homogeneous, we decompose the perturbations into waves as described in the previous 
section. The only wave that can propagate upward is the fast magnetosonic wave /— . The boundary condition requires 
that its amplitude is zero: 

^V" = 0. (A37) 

Computation of the cycle efficiencies 

The coupling efficiencies at the shock are computed by decomposing into waves the perturbations just below the 
shock. For a given value of u), these perturbations are set by the upper boundary condition described above. The 
coupling efficiency Q l sh is then computed as the ratio of the amplitude of the wave i (i — s±, f+, a±, S) to the amplitude 
of the upward propagating fast magnetosonic wave: 



SA 



(A38) 



sh 



To obtain the coupling efficiency Qy, we need to determine in what proportion an upward propagating fast wave 
(/— ) must be added to a wave i at the upper boundary at z s h, so that the lower boundary conditio n is r espected. For 
this purpose, we successively set the perturbation at the shock to : (i) a wave i (described in table IA.4j) of amplitude 
SA l sh = 1, (ii) an upward propagating acoustic wave (/— in table PO| of amplitude 6AL~ = 1. We then integrate the 
differential system from the shock to z\ ow = z\> — 3iJy j where the lower boundary condition is estimated through the 

amplitude of the upward propagating wave /— wave: x = (j^J ■ This resulting quantity is called x 1 for the wave 
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for the case (i) and x* for case (ii). The leaking boundary at z\ ow is respected for a combination of the waves i and 
/— (at the shock) if : 



5A^xf-+5Al h x* = 0. 



(A39) 



Thus the coupling efficiency can be computed as : 



n i _ $A l sh _ x l 



(A40) 



VERTICAL MAGNETIC FIELD 



In this appendix we write the equations governing the perturbations in the case of a vertical magnetic field, using 
the same method as described above in the case of a horizontal magnetic field. 

Definition of the new variables 



Sh 

Sf 
SE X 

Sv x 



5R 



v 



Sp _|_ Svz 

p V 

vSv z + Mf 
6v x - v 5 -%- 

x v B 
c SB V 

SVy - 



SVy = SVy 

8S 



v B 



7-1 V P 1 P ) 



Sp 
■ 

Sv z 

V 

SB Z 
B ~~ 

Sv x = 

B 

SVy = 
SBy _ 

B 



1 

1-M 2 
1 

1-M' 2 



^-M 2 5h-SS 

- 5 -l + SEy + 5S 
^SE X 

(Sv x - 8E X ) 



^5] 

Differential system 



(SVy ~ SEy) 



dSh 



dz 



OSEy 

dz 

dSvy 
dz 

dSS 

dz 



dz 


V 


(1-M 2 ) 


dSf 




iujv 




dz 


1 


-M 2 




dSE x 




iuj 


dz 


V 




dSv x 







% - M 2 Sh - SS 



% + 5h+ l - - ; + 



i \ ss' 

7 . 



5v x - ^5E X 



M 



V (l 



[SE X - Sv x ] , 



SE X 



+ 



ik x v 
1-M 2 
iuj 



iu> 



-k + 5h + (7-1 + T72 ) — 



M< 



7 



[SEy - SVy] 



(Bl) 

(B2) 
(B3) 



(B4) 
(B5) 

(B6) 

(B7) 
(B8) 
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TABLE 2 

Decomposition into waves: vertical magnetic field 





Fast ± 


Slow ± 


Aliven ± 


Entropy 


Sh 
Sf 
5E X 

Sv x 

SEy 
SVy 

SS 


c 2 A + k+v\ 







Sp 
P 

Sp 
P 

, k z v \ Sp 
~ J ~J 

V A k z v \ Sp 

~V T Oc J P 


(l + k£L) Sv x 












(l±^)<Ha± 



(l-7)f 
C 2SS 
1 







ss 


k z 


to M^f-p _j_ »| k t c 
c 1 — fA 2 2c' 2 fxui 


to _j_ 1 V A to 3 u 
^^i;a 2 v '^ c 2 k 2 v 2 +L*j 2 




V 



Wave decomposition in a uniform flow 
Boundary conditions at the shock 



Sh 



sh : 



Sv xsh ■ 

SVy S h - 
1 



iuA( 

V 2 



£/ sh = iawiAC 1 



V2 
Vl 



Sfsh = ik x V!A( 1 



= 0. 

k x 

= o, 



V2 
Vl 



V-2 
Vl 



d<P AC 



£2 
Vl 



(B9) 

(BIO) 

(Bll) 

(B12) 

(B13) 
(B14) 

(B15) 
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